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Abstract 

A detailed investigation of the theoretical ambiguities present in the QCD de- 
scription of photon production in e + e~ annihilation is given. It is pointed out that 
in a well-defined perturbative analysis it is necessary to subtract the quark-photon 
collinear singularities. This subtraction requires the introduction of an unphysical 
parameter in the perturbative part of the cross section. The subtracted term is 
factored into non-perturbative fragmentation function. The dependence on the un- 
physical parameter cancels in the sum of non-perturbative and perturbative parts. 
It is pointed out that for E y < \/i/(2(l + e c )) the non-perturbative contributions are 
suppressed. Using a general purpose next-to-leading order Monte Carlo program, 
we calculate various physical quantities that were measured in LEP experiments 
recently. 
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1 Introduction 



The production of a photon (or an isolated photon) in association with hadrons in e + e~ 
annihilation is a useful process to learn about the differences in the properties of q(pi 
and qqg final states, to measure the parton-photon fragmentation function and to test 
QCD predictions in a channel crossed to photon-photon annihilation. The corresponding 
theoretical problems are well understood in the case of prompt photon production at 
hadron colliders, photo-production of jets and heavy flavor and photon-photon scattering. 
It is an important development that experiments at LEP give us high statistics data and 
open ground to study even photon plus multijet final states [l], The better data call for 
a quantitative QCD description. 

The QCD description of inclusive photon production has a simple, but important fea- 
ture: the photon has hadronic component. In the perturbative treament this fact is 
reflected by the appearance of collinear photon-quark singularities. In order to obtain 
well defined cross sections in perturbative QCD in all orders of the running coupling a s , 
these singularities are to be subtracted and absorbed into the photon fragmentation func- 
tions (factorization theorem) || ||. The fragmentation functions of the photon satisfy 
inhomogeneous evolution equation; it grows with Q 2 therefore, it is called "anomalous" 

It is also interesting to study the case of isolated photon. Physical isolation means 
that we isolate the photons from hadrons and so we cannot make distinction between 
quarks and gluons. Gluons, however, cannot be isolated completely from the photon 
without destroying the cancellation of soft gluon singularities between the virtual and real 
gluon corrections. Therefore, a physical isolation cannot eliminate completely the collinear 
photon-quark singularities, and so, even in the case of isolated photon production the cross 
section contains "anomalous" (non-perturbative) piece. This problem has been recognized 
clearly in the next-to-leading order QCD study of isolated photon production at hadron 
colliders || |[. The theoretical subtleties of defining isolated photon cross section in 
perturbative QCD, however, have not been clearly formulated in previous studies in the 



case of e + e annihilation [1C, 11 



In section 2 we review the next-to-leading order description of the inclusive (non- 
isolated) photon production. In section 3 we outline the change in the formalism due to 
the introduction of isolation cuts for the photon production. We point out that isolation 
cannot completely eliminate the non-perturbative fragmentation contribution, although it 
can reduce its size. In section 4 a detailed perturbative study is given for the cross section 
of isolated photon plus jet production up to order 0(aa s ). We review the mechanisms 
of the cancellation of the infrared singularities and point out that in perturbation theory 
for processes containing a photon in the final state the definition of a finite hard scatter- 
ing cross requires a counter term which necessarily introduces an unphysical parameter. 
Section 5 contains our numerical results for the isolated photon plus n-jet production at 
LEP. To demonstrate the flexibility of our numerical program to calculate any jet shape 
parameters, we calculate the distribution of the photon transverse momentum with respect 



1 



to the thrust well. The last section contains our conclusions. 



2 Inclusive photon production in e + e annihilation 

According to the factorization theorem, the physical cross section of inclusive photon 
production is obtained by folding the fragmentation functions D 7 / a (x, [if) with the finite 
hard-scattering cross sections da a : 

da ^ rVs/2 ri d<7 a 

-7-^ = 22 dE <> / dxD l/a (x^ } )—^{E ai ^^f 1 a s {^))5{E 1 -xE a ) 1 (1) 

a-c/ 7 a Jo Jo ci-c/ a 

where a s (/x) is the strong coupling constant at the ultraviolet renormalization scale [i and 
\if is the factorization scale. 

It is instructive to investigate the decomposition of this generally valid expression up 
to next-to-leading order. First we remark that 

D lh {x) = 6(1 -x) + 0(a 2 ), (2) 

therefore, to leading order in the electromagnetic coupling, the term in eq. ([!]) given by 
a = 7 is a purely perturbative contribution. We use this equation to eliminate D 1 / 1 (x) 
from eq. ([!]). The hard scattering cross section da^jdE^ is of order a in comparison to the 
leading order annihilation cross section ao-Q The leading non-perturbative part given by 
the fragmentation function, however, is of order a/a s . This contribution is the "anoma- 
lous" photon component. Its enhanced order is due to the fact that the scale dependence 
of the fragmentation functions D 7 / a (x, ///), a = q,q,g are given by the inhomogeneous 
renormalization group equations [12], [I 



V^- D j/a(x, fj) = -P-y/a(x) + — ^ / — D i/b I — , A* 1 Pb/a(y), (3) 

where P&/ a (x) denote the Altarelli-Parisi splitting functions. To order aa s the inhomoge- 
neous terms have the expressions {§] 

Pl/a(x)=P!}%{x) + ^f a (x) 1 (4) 



where 

X 

and after trivial replacement of the color factors in eq. (12) of ref. 0, we have 



P» W ^i±(lzf)! p<o. (l) = , (5) 



pS M = (6) 



e 2 CV \ —- + -^x + ( — 8 + -x ) log x + 2x log (1 — x) + ( 1 — -x ) log 2 x 











r + (- 8 


+ 5*) 









+ 



4 

log 2 (1 - x) + 4 log x log (1 - x) + 8Li 2 (l - x) ix 2 
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1 In the following analysis, when the order of a contribution is given, it is always understood in 
comparison to the leading order annihilation cross section ctq- 
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P?fa) = (e 2 q )T R {-A + 12x-^ + fx^ 

16 16 \ 1 
10 + 14s + — x 2 + -x- 1 ) log x + 2(1 + x) log 2 x \ . 



In the last equation, 

Np 

<<3> = £4 (7) 

9=1 

where iV^ is the number of flavors. 2 

The unique solution of these inhomogeneous equations requires non-perturbative input^] 
at a certain initial scale //. At asymptotically large values of /i, however, the solutions are 
independent of the initial values and one obtains 

(X /i 2 

Q/ /i 2 

J^^/sO^) = 2^ lo ^-^ a i/g{ x )- ( 9 ) 

Exact analytic expressions for the Mellin transforms of the a a / 7 (x) functions have been 
found in refs. [^|, |5[]. These are related to the a 7 / a functions via crossing. It is useful, 
however, to have a parametrization in x-space. Formulas which accurately reproduce the 
exact leading logarithmic solutions were given in ref. 



<h/*( X ) = e l\ 



2.21 1.28x + 1.29x 2 049 nnon „.\2„.-i.54 



1.631og(l -x) 



-x au4y + 0.002(1 - x) V 



(10) 



a 7/ff (x) = ^[0.0243(1 -x) L03 x- a97 ]. (11) 



A new parametrization of the photon fragmentation functions is described in ref. ]20| . The 
most striking feature of these solution is that they increase as l/a s with increasing the 
evolution scale. Therefore, at high energy the contribution from the quark fragmentation 
into a photon gives the leading order (a/a s ) term 

do-(°) 4 e 2 (IE \ 

dir^?<^'«&") +0(o) ' (12) 

In next-to-leading order, the \i dependence of -D 7 / 9 has to be calculated with the next-to- 
leading order evolution equation and we should also add the order a hard scattering cross 
section 

d<7 7 



d£ 7 




2 We assume e + e~ annihilation via virtual photon. In order to obtain formulas valid at the Z° peak, 
trivial modifications of charge factors are required. 

3 In the literature it is usually called Vector Meson Dominance (VMD) contribution |L3[ |l4| || . 
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+ V / dE a / dxD^/ a (x,iif)-r^ r (E a ,iJ l ,iif,a s (ii))S(E J -xE a 
a Jo Jo a.c/ a 



+ -r^-(-E' 7 ,/i/,a s (/i)) + 0(aa s 

Q-C/-V 



(13) 



where da^ /dE a denotes the order a s cross section of quark and gluon production. 

The 0(a) hard-scattering cross section da^ /dE^ is defined by subtracting the photon- 
quark collinear singularity in the MS scheme 



djff 



lim , 



'dcrW d<7$ 



+ 



CT 



d£L 



(14) 



where the first term on the right hand side is the partonic cross section in 4 — 2e dimensions 
as defined by Feynman diagrams 



dcr( > _ v e\ a 2 



'47T/X 



2\ e 



1 



Hi-el 



x 



V2/13 2/23/ 2/132/23 



dyi2 dyi3 d?/23 0(1 - 2/i3 - 2/23) (2/122/132/23) E 

(15) 

2EL\ 



5(1-2/12-2/13-2/23)^^1-2/12- /s 



where H = 1 + while the second term is the MS counter-term 



da, 



(0) 
CT 



(LEL 



27r£r(l -e) v^ ' /l! 



dM°) 



^ v ; d£„ 



■{E q )8{E 1 - xE q 



(16) 



The integrations in eqs. fllBD , fllq) are easily performed. The collinear poles cancel in their 
sum. Setting e — 0, one obtains 



daf 



o 7= L 7i\ P 7/g x 7 lo § 2 



(17) 



where x 7 = 2E 1 /\js. 



The C(a s ) corrections to the da qj9 hard-scattering cross sections are defined by the 
Feynman diagrams of fig. 1. First we note that da^ can be obtained from da^ by 
modifying the charge factors: 



daW 
dE n 



a x 4 



'S(l - Xg)X\ 



It 



The cross sections da^ receives both real and virtual corrections. The loop correction can 
be written as 



daioop „ oi s 2 1 n 

— — = C F a - -=ti - 

dE q 2tt y/S \ s 



2\ £ 



(4tt) 



em-e) 



3 + (vr 2 - 8)e 



5 



E n 



(19) 



The Bremsstrahlung contribution has an expression similar to da^ /dE 1 (|IJ): 



der, 



real 



dE n 



x 



C F a 



e 2 q a s 



{el) 2vr ^ 



H 



2\ e 



1 



2/23 3/13 \ 
,2/13 2/23/ 



s J r(l-e 

%12 - £2/132/23 
2/132/23 



d2/i2 dyi 3 d?/23 0(1 - 2/i3 - 2/23) (2/122/132/23) 



(20) 



5(1 - 2/12 - 2/13 - 2/23)5 1 - 2/23 - 



2E 1 



The sum of the loop and Bremsstrahlung contributions has the simple expression 



do« 
dE n 



C F a 



(47T) 



2^^ er(l-e) 



(21) 



x 



'log(l 



1 Xn 



3 
2 



x 



1/ + 



2 

3 5 

-X a H 

2 q 2 



l + xl 
2\ogx q - q - 



where the index + denotes the usual "+ prescription" of regularizing singular behavior at 
x q = 1. The remaining single pole is cancelled when one adds the MS counterterm da^ T 



which is defined as 



dcr, 



(i) 

CT 



dE n 



C F a 



(4tt) 



P 



(o) 



(2E n 



2tt ^seT(l-e) q/q \^ 



(22) 



The final result is obtained after setting e = 0: 



da (1) a 2 

W^o- 7= 

27T VS 



d£„ 



(23) 



:i+x 



log(l 



+ 



X n 



-7T 



5(1 
1 



X q ) + 2hgX q 



l + x. 



Xn 



l-X n 



2 V 1 - x a 



2 Xq+ 2 



where x q = 2E q /y/s. This result can also be deduced after replacing trivial color factors 
from the coefficient functions of inclusive single hadron production calculated in ref. ||. 

The theoretical input described in this section is sufficient to extract the photon frag- 
mentation functions from experimental data in next-to-leading order accuracy. A complete 
analysis requires the measurement of the inclusive photon production cross section at var- 
ious energies. The recent LEP data give information at the Z-pole. Unfortunately, the 
data obtained at PETRA, LEP and TRISTAN suffer from low statistics. Needless to say 
that such an experimental study would give very important complementary information 
on the fragmentation functions of the photon obtained at hadron colliders. 
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3 Inclusive isolated photon production in e + e anni- 
hilation 

Let us now consider the inclusive photon cross section with photon isolation. One can argue 
that due to isolation cuts the fragmentation contribution is suppressed. As a consequence, 
isolation changes the relative importance of the different contributions. It is reasonable to 
consider the effect of isolation typically as an order a s effect. After imposing the isolation 
cuts, the fragmentation contribution will be of order a, i. e., the same order as the order of 
the pointlike perturbative cross section daj^/dE^. Isolation in practice can only be made 
with finite energy resolution. Therefore, we require that in a cone of half angle S c around 
the photon three momentum the deposited energy be less than a fraction e c of the photon 
energy. In experiments this parameter e c has a value typically about 0.1. Calculating 
da^ iso /dE 7 , we should insert a combination of 6 functions in the phase space integrals as 
follows 

S(e c ,5 c ) = 9($ qi -5 c )9($ qi -5 c ) 

+ 9($ qi -5 c )9(5 c -$ qi )9(e c E 1 -E q ) (24) 
+ 9($ qi -8 c )9(8 c -$ qi )9(c c E 1 -E q ). 



Let us require that 
and choose integration variables 



1 , . 2 5c 1 

e c < - and sin — < - 

2 2 2 



_2E 1 _ yia 

Xry , y 

x 7 



We define the hard scattering cross section again with a collinear counter-term 

^ = lim(^ + ^^V (25) 
da; 7 e^o I dx 7 dx 7 J 

where the first term in the right hand side is calculated as given by Feynman diagrams in 
4 — 2e dimensions and for the counter-term, we use the MS-type expression 

d4°T,iso _ a v e\ (Any (0) 1 

" 2(T °27 \ yi)eT(l-e) P ^ )0[X " ~ 1+1/ ' (26) 

After performing the integration over y and setting e — 0, one obtains 



d<ri°L _ a ^ e 2 n fr„rm, _ s(l - x^xiy 



' 7. is" _ o a _ ~ q 



9 x 



l + e c/ 
(27) 
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+ P$(* 7 )log - e 2 qXj (l - 2y c ) 



where y c and y m are defined as follows 

1 — x 7 . 2 • J i i 1 

1 - x 7 sin ^2 [ x n 

One can make several comments on this result. 

• The unisolated case can be recovered in the limit e c — > oo (cf. eq. (|T 

• Imperfect isolation allows for a contribution from the fragmentation: the photon 
looks isolated since the relatively soft fragments surrounding it are not counted. 

• Assuming perfect energy resolution (e c = 0) we obtain vanishing counter term. In 
higher order, however, we can not isolate the photon from the soft gluons completely 
(we shall discuss this point in great detail in the next section), therefore, one can 
not set the value of e c to zero. 

• In the leading logarithmic approximation one can define a fragmentation function 
with isolation satisfying a modified inhomogeneous evolution equation 

d 

fJ-Q^D^/ a (x,fj,,e c ) = (28) 

a , ( 1 \ a s f dy fx \ 

-P l/a {x)6 x - — — + — 2^ / — D l/b -, [i, e c P b /a{y)- 

Clearly, if D 1 / a (x,n) is a solution of the evolution equation without isolation then 

D y/a (x, /I, e c ) = D l/a (x, fi)6 (x - 1 + g ) (29) 

will be the solution of the evolution equation with isolation. In next-to-leading 
logarithmic approximation and/or choosing a different counter-term (for example 
completely subtracting the contribution of the singular region as defined by the third 
term of eq. (|23|)), the isolated fragmentation can also be dependent on 5 C therefore, in 
general, one cannot simply identify the isolated fragmentation with the non-isolated 
fragmentation in the high-x region. 

In next-to-leading order, the physical cross section of isolated photon production is given 
by the terms as follows 

d <V x n _ , d4o } , -a 2E, 2E, 1 

d^ (ec ' 5c) ~ d^7 + d^7 + 2(70 2^ ( 7^ ' ^ ( 7i " tttJ 

+ £/ dE « / 1 ^D^ a {x,yLf)j^{E ait i, l i f ,a a iji))S{E ri -xE a ). (30) 
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This decomposition is scheme dependent. The first term on the right hand side of this 
equation has been calculated in the MS scheme (see eq. (p7|). It also appears useful to 
calculate the next-to- leading order perturbative cross section da^ /dE^ in the MS scheme. 
This requires the calculation of the next-to-leading order splitting function P^) a in the 
presence of isolation cuts and a corresponding local subtraction term has to be found. 
This is a complex but feasible calculation. Since such a result is not yet available, in 
the next section we carry out the calculation of da^ in a different subtraction scheme 
where the photon is completely isolated from the quarks but not from soft gluons ( "cone 
subtraction"). In this scheme, in leading order, the counter-term is vanishing and the cross 
section becomes independent of e c : 

d<r<°> „ a ^ el f (0 , 1-y. 



^ - ^ T, ^ [F*?.M in -*,<!- 2 V ,) J . (3!) 

We note that the logarithmic divergence at x 7 = 1 is the usual soft singularity. Contrary to 
the case of the MS scheme, with cone subtraction the cross section is continuous and always 
positive (see fig. 2 for comparison). One may argue that in this scheme the perturbative 
part is separated more efficiently, consequently the contributions of the non-perturbative 
terms (proportional to D™J a ) become relatively smaller. 

In general we find that the non-perturbative terms contribute mainly in the region 
x 7 > 1/(1 + e c ) thus we conclude that the perturbative predictions appear to be reliable 
for £ 7 < v/i/(2(l + e c )). 

In the next section we present the results of our next-to-leading order perturbative 
calculation of da^o f° r isolated photon plus n-jet production. We conjecture that a jet 



algorithm applied to the isolated photon hard scattering cross section (eq. (|30|) ) provides 
an infrared safe isolated photon plus n-jet cross section. This is supported by the fact that 
our isolation prescription does not influence the soft-gluon structure of the cross section. 
If we can define a jet algorithm, 

j-iso Hrr iso Hrr iso Hrr iso 

^ *> = ^ «» + + %f «> + ' ' ' ■ < 32 > 

such that every term on the right hand side is finite and we count every particle only once, 
then isolated photon plus n-jet cross section appears to be infrared safe. 

We shall see that the non-perturbative ( "anomalous" ) contributions are important only 
in the case of photon plus 1 jet production when the cross section is dominated by the 
x 7 > 1/(1 + e c ) region. 

4 Isolated photon plus n-jet production 

In QCD, the differential cross section at 0(a 2 s ) is a sum of the real and virtual corrections: 

da = |M 4 | 2 d^ (4) + |M 3 | 2 d^ (3) , (33) 
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where dS^ is the n-body phase space element with the flux factor included. For infrared 
safe quantities both terms on the right hand side are separately divergent, but the sum 
is finite. It is very difficult to handle numerically this cancellation. Fortunately, at least 
at one loop, the divergencies can be cancelled analytically. There are two commonly used 
algorithms to achieve such a cancellation — the subtraction method [|16| , [17| and the phase 



space slicing method JL8], [L9|]. They both rely on the fact that after partial decomposition 
|M 4 | 2 can be written as a sum of terms with single pole singularity. Focusing our attention 
to the case of qcpig final state, we find four such terms: 

W = c F aa, (Mw + M?l + M™ + Mll\ (34) 
V Vgq Vgg Viq Vig ) 

where 

Mi = (Pi + Pjf/s, = M|) (35) 
The pole part of each term is defined as 

P 

where P i7 - = lim M ti . (36) 

It can be integrated analytically over either the whole or a part of the phase space. In this 
way, in general, we obtain analytical expressions for the regularized divergencies of da 1 - 4 -* 
which cancel against the divergencies of the virtual corrections, da^ (KLN theorem). 
When a photon in the final state is observed, the cancellation mechanism described above 
does not apply to the y iq (q) poles. The reason for this is that the process is exclusive in 
the photon and the virtual corrections with the photon in the loop cannot contribute for 
kinematical reasons. 

To make the discussion more transparent, let us consider contributions from the region 
where only y qg is small. The virtual corrections can also be split into three terms 

|M 3 | 2 = C F aa s (Mj + Mgf + m} 3) ) , (37) 

such that contains one half of the singularities, the second term contains the other 

half and the third is the finite part.0 (Notice that there are no M$, M§j terms.) Then 
we shall concentrate on 

Mllas^+M^dS^ (38) 
Vij 

parts of the cross section. 

In the subtraction method, one considers the combination 

Mil dS W - SudS^Q + ( I S^dS®) dS^ + M®>dS®, (39) 
Vgg Vgg V Vgg J 



4 For the reader's convenience, we give the explicit expressions for My, and in the appendix. 
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where the integration over dS^ is meant to be an integral over the gluon variables. dS"^ 3 -' 
means the factorized four-body phase space element in the limit when the gluon is soft or 
collinear to the quark: dS^^ = dS^dS^. 

In the phase space slicing method, formula (BSJ) is written as 



^6(y gq - y )dS^ + - y gq )dS^ + M^fdS^. (40) 

Vgq Vgq 

If y is chosen small enough (y < 10~ 4 ), then 

^e(y gq - y )dS^ + ^9(y - y gq )dS^ + M$dS® (41) 

Vgq Ugq 

is a good approximation. The first and second terms depend on y strongly, but their sum 
is independent of this unphysical parameter. The strong y dependence originates mainly 
from the slicing of the soft gluon region. 

If one wishes to calculate isolated photon production, one has to make sure that the 
restriction of the phase space does not disturb the cancellation mechanism of soft and 
collinear gluons. At hadron level the meaning of photon isolation is well-defined. At 
parton level however, one has to be careful because the isolation prescription is different 
for states with different number of partons.Q If the photon is isolated form the partons 
with y 7 , we should include isolation cuts with respect to all partons: 

%7? - VjWivjg - 2/7)%79 ~ V~f)- ( 42 ) 



The isolation from the gluon can be implemented only in the first term of formula (£ 
However, if we cut the soft gluons in the first term, then the cancellation of singularities 
between the first and second terms breaks down. One can maintain the cancellation 
introducing an energy resolution parameter e such that a gluon is isolated from the photon 
only if its energy is greater then eE 1 . Accordingly, isolation for the first term means 
multiplication with 

% 7 g - y~/W(y~/q - y~i){ 1 - % 7 - y~/g) e ( E 9 - eE v))- ( 43 ) 

Clearly, this criterium is "not physical" in the sense that one cannot implement it at 
hadron level since we apply different cuts to quarks and gluons. 

If we introduce photon isolation in the slicing method, from formula (f41~D we obtain 



— 0{y gq - y )dS^e(y, q - y,)6(y iq - y y )9(y jg - y y )+ (44) 
Vgq 



5 One may object isolation at parton level arguing that the fragmentation process inevitably scatters 
hadronic matter into the isolation cone. For a purely perturbative analysis, this objection is not valid. 
To understand the reason for this, let us consider the same process at higher energies, say ^fs = lOTeV, 
in which energy region perturbation theory is expected to give even better description. Clearly, at this 
energy, fragmentation does not alter the energy flow, therefore isolation at parton level corresponds to 
isolation at hadron level. 
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(^6(y - y gq )dS^ + M®dS<A 9(y 7q - y 7 )9(y 7q - - y 7 ). 
\ Vgq / 

Usually, y 7 ^> y . This means that when changing y at fixed y 7 , the contribution from the 
soft gluons will be cut independently of y and consequently, the y dependence is damped 
in the first term. On the other hand, in the second term the yo dependence is not damped 
by the gluon-photon isolation condition. The conclusion is that the yo dependence will be 
different in the two terms and, therefore, the cross section will depend on the unphysical 
parameter y . It is important to notice that if y gq > yo, then there exist an e' such that if 
e < e' then 

1 - 9{y 7 - y 7g )6(E g - eE 7 ) = Q{y ig - y 7 ), (45) 

therefore (f|4|) defines a finite cross section, but y plays in a sense the role of e used in 
formula flUjp . 

To demonstrate the y dependence of the isolated photon cross section explicitly, we 
calculated the isolated photon plus 1- and 2-jet cross sections using the isolation criterium 
given by formula ([44]). First we make two technical remarks about the slicing method. 

Choosing large yo, the pole approximation is not precise enough in the singular region; 
one has to take into account the non-singular terms in the same region, i. e., one should 
add the 

(MmasM - ^EldS^A % - Vgq) (46) 

V Vgq V 9 q J 

correction term. The calculation becomes analogous to the subtraction method and one 
has to introduce the e energy resolution parameter. 

It is a practical question to establish at what value of yo the correction term fl46|) 
becomes important. The most straightforward way to calculate the finite terms in ( |14] ) is 
to perform the integration by a Monte Carlo method, which leaves sufficient flexibility to 
calculate any jet shape parameter one wishes to obtain. The Monte Carlo calculation has a 
finite statistical error which of course, can be reduced by generating more points. Then the 
criterium which determines the importance of the correction term (f46"D is to require that 
the systematic error introduced by neglecting (|4*6| ) has to be smaller then the statistical 
one. Clearly, this critical value of yo depends on the jet resolution parameter yj as well 
as on y 7 . For the case of 3-jet production without photon in the final state, an analysis 
was carried out in ref. [[0| to determine the critical value of yo above which the systematic 
error dominates. They found that choosing yo/yj < 0.01 removes the systematic error. 

The number of isolated photon plus ra-jet events can be conveniently parametrized in 
the form 

— a 7+niets (yj, y ) = ^Y,T^\ Vo) + TT-gSHvJ, Vo)) (47) 

4 

= ir H -rk9n(yj, i/o)(i + a s R n (yj, yo)), 

Z7V q \ e q) 
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where o"o is the leading order cross section of the reaction e + e~ — > hadrons and the 
Rn(yj,yo) functions are defined by the equation. In figs. 3 and 4, we show the yo de- 
pendence of the 0(a s ) QCD corrections, R n (yj,y ), to the isolated photon plus 1-jet 
and the isolated photon plus 2-jet cross sections. To obtain the corrections, we used the 
following algorithm: 

1. select isolated 7 + n-jet events by requiring the invariant mass of the photon with 
any particle in the event to be larger than y 7 (see formula 



2. apply EO cluster algorithm to the hadronic part of the event; 

3. separate 7+ 1-, 2-, and 3-jet events by the number of remaining clusters of hadrons. 

We used y 1 = yj. As expected, the y dependence in R n {yj,yo) is strong up to y = y 1 . 
As explained before, in formula ([44]) the yo cut plays the role of the e parameter of formula 
(fj3"|). Therefore, the (apparently) physical cut, ( fP| ) is in fact unphysical because yo is no 
longer a dummy variable of the cross section. 

In order to demonstrate that we control the numerical evaluation of the integrals at 
small yo values, we calculated R n (yj,yo) with 6{y ig — y 7 ) in ((PJ) removed. We denote the 



corresponding quantity with R n (yj,yo). According to the discussion after formula 
this alteration should remove the yo dependence. The explicit calculation shows that this 
indeed happens.^ We see that in order to obtain a yo independent result, we have to use 
an unphysical cut: different cuts are applied to quarks and gluons. 

We conclude from this discussion that if we want to define a finite isolated photon 
plus n-jet cross section we have to make a subtraction which depends on some unphysical 
parameter no matter which algorithm — the subtraction or slicing one — is used (see 
formulas ([43"[), ([44]) and ([45])). In other words, the physical isolated photon plus n-jet 
cross section always contains some non-perturbative ("anomalous") contribution which 
is expected to give contributions comparable or somewhat smaller than the 0(a s ) QCD 
corrections. For the separation of perturbative and non-perturbative parts of the cross 
section, one must introduce an unphysical (non-zero) parameter. Of course, the sum of 
the perturbative and non-perturbative pieces is independent of this parameter. In the 
previous section we pointed out that the non-perturbative contribution is expected to be 
small for x 1 < 1/(1 + e c ). 



5 Numerical results 

As advocated in section 4, we carry out our calculation with the subtraction method. The 
event definition is the following: 

6 In fact, we can see weak y dependence in R n (yj,y ). The origin of this dependence is the use of the 
pole approximation. It can be observed for values yo > 10 -3 (this value depends on y,j) in accordance 
with the observation made in ref. ]l9| . The results are shown for yj = 0.06. For other values of yj the 
dependence is similar. 
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1. isolate the photon; 



2. apply EO cluster algorithm to the hadronic part of the event; 

3. separate 7+ 1-, 2-, and 3-jet events by the number of remaining clusters of hadrons. 

Photon isolation can be achieved either by isolating the photon in a cone (cone isolation) 
or by requiring the invariant mass of the photon with any particle in the event to be 
larger then an invariant mass cut y 1 . From experimental point of view the cone isolation 
is more natural. Unfortunately, the results by OPAL are corrected experimental values 
in order to compare the measured rates with the matrix element calculation of [|10j where 
invariant mass isolation was used (with y 7 = yj). We give results for both. Since the QCD 
corrections are very sensitive to the event definition we give explicitly how formula (|3U|) is 
modified in the case of cone isolation: 

e{$ qi - 5 c )e{$ qi - 5 C ) 

x \(l-9(5 c -# gj )9(E g -e c E 1 ))^-dSW (48) 
I Vgq 

_^El dS (^3) + ( f £il dS ( 9 )\ dS (3) + M (3)l 

Vm V Vm J J 

and in the case of invariant mass isolation: 



% 97 - yj) d (y^ - vj) 

x !(l-e(yj-y g7 )9(E g -e c E^))^dS^ (49) 
I Vgq 

L d5 (4-3) + ( r En dS ( 9 )) dS m + M (3)l 

Vgq J J 



Pgq. 



Vgg 

To obtain the isolated photon plus n-jet rates, the formulas above are multiplied with 
9 functions as follows: 

• One photon plus 3-jet: 

0(y qg - yj)9(y qg - yj)9(y qq - yj)9(y qi - yj)9(y qi - yj)9(y gi - y.j). (50) 



• One photon plus 2-jet: 

Denote i and j the partons which when combined have the smallest invariant mass in 
the hadronic part of the event, so they are combined into pseudoparticle c. Denote 
k the third parton. In the three-body phase space the momentum of the j particle 
is identically zero. Then we use 

%j - VijWiVck - yj)0(y n - yj)0(y kl - y.j)- (51) 
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• One photon plus 1-jet: 

%j - y qg )0{yj - yq 9 )0{yj - y q q)0{y.j - y<*)- (52) 

In the case of cone isolation, we also required that the energy of the photon has to be 
larger than 7.5 GeV. The half-opening angle of the cone is 15°. 

We shall give the results of our calculation for the partial widths T(Z — > 7 + n jets) as 
ratios to the hadronic width: 

r(Z^ 7 + Ilj ets)_ + Mi m 



T Z 



hadrons) " ^ + 3q) (1 + Si + L42 (^.) 



where 



c f = v) + a} (54) 
and Vf and aj are the weak vector and axial vector couplings: 

v f = 2I 3J - Ae f sm 2 9 w (55) 

a, = 24/, (56) 

so with sin 2 9w = 0.23, v u = 0.39, = —0.69, a u = +1 and ad — —1. The g n (yj) functions 
can be expanded in a s : 

9n(yj) = $\yj) + ^9 { n\yj) = + oc s R n { yj )\ (57) 

and our aim is to compute the g$ functions (of course, g^\yj) = 0.) 

It is interesting to study the photon energy spectrum of the jet cross sections 

j iso j iso 

^g^c,e c ) (58) 

at some realistic values of the isolation parameters 5 C , e c . In fig. 5 the photon energy 
distibutions of one jet and two jet production are shown in the Born approximation, while 
in fig. 6 the same curves are plotted but including the next-to-leading order corrections. 
Due to obvious kinematical reasons one jet production is completely dominated by the hard 
photon region x 7 > 1/(1 + e c ). In this region as we pointed out one may get substantial 
(but not overwhelmingly large) "anomalous" photon contribution. It is difficult to estimate 
the "anomalous" contribution since we do not have yet enough phenomenological input. 
Certainly a combined study of the hadron collider and LEP data would help to understand 
its size better. We note that in the high x region the application of perturbative QCD by 
itself requires some care due to the appearance of large logarithms of type log(l— x). Indeed 
the QCD corrections are larger for one jet than for two jet production. It is interesting 
to compare the one jet data to the perturbative QCD prediction, but one should not 
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be surprised if one does not find perfect agreement. The non-perturbative corrections 
appear, however, negligible in the case of 2-jet production since it is dominated by the 
complementary region x 7 < 1/(1 + e c ). Requiring that < 1/(1 + e c ), the 2-jet results 
remain practically unaffected, while this cut largely eliminates the 1-jet production. This 
is illustrated by the numbers given in Table 1. There is a tendency that if we shrink the 
isolation region the perturbative contribution increases. From figs. 5, 6 and, we can see 
also that the total one jet and two jet rates should depend weekly on e c . The reason for 
this is that the photon energy distribution changes weakly if we change e c in the physically 
interesting region of 0.06-0.2 . 

In addition to the ambiguities due to "anomalous" photon production there are also the 
usual scale ambiguities. In fig. 7 we present the predicted values of the T(Z — > 7 + njets) 
(n = 1,2) partial widths for the cone isolation with e c = 0.1. The bands between the 
dashed lines represent the scale dependence between the scales M z /2 and 2M Z - We used 
a s (Mz) = 0.12 and a = 1/137. The e c dependence is so weak for experimentally feasible 
values that the uncertainty introduced by the e c dependendce is much smaller than the 
scale dependence and therefore we did not show it. The scale dependence of the 1-jet rate 
is rather large. This is a reflection of the fact that the QCD corrections are large. In figs. 
8 and 9 the same curves as in fig. 7 are depicted in the case of invariant mass isolation 
with y 1 = yj for the 1-jet and 2-jet rates, respectively. In the same figures, we show the 
enhancement induced by the choice of smaller isolation region. In accordance with our 
previous discussion, the enhancement is larger for the 1-jet rate than for the 2-jet rate. 
We note, however, that when comparison is made to the data at a given isolation it is 
important to use exactly the same isolation and event definition both in the experimental 
and theoretical analysis. Therefore one can not just change the value of y 1 such that the 
prediction fits better the data. In particular one is not allowed to use different values of 
y 1 in case of one jet and two jet production. In a given subtraction scheme with well 
defined experimental isolation cuts all the parameters of the perturbative part are fixed. 
In particular the discrepancy between the measured 7+ 1-jet rate and the perturbative 
prediction at y 7 = yj may indicate non-negligible anomalous contributions. 

As mentioned in the section 4, the Monte Carlo approach is useful because it leaves 
sufficient flexibility to calculate any jet shape parameter. To demonstrate this feature of 
our work we present the result of matrix element calculation for the distribution of the 
photon transverse momentum with respect to the thrust axis (fig. 10). The thrust axis has 
been calculated all particles taken into account, including the photon. We used invariant 
mass isolation (with y 1 = 0.005 and 0.06) to isolate the photon from the partons. We 
also required the photon to be more energetic than 7.5 GeV. For small p?, configurations 
with thrust value close to one may occur. The histogram is normalized to one, therefore 
the uncertainty in the small pr region influences the behaviour in the large pr region. We 
note, however, that requiring x 1 < 1/(1 + e c ) the small pr region will be suppressed. 

Finally, in fig. 11, we present the predicted values of the Y{Z — > 7 + n jets) (n = 1, 2) 
partial widths for the cone isolation with e c = 0.1 when Durham clustering algorithm is 
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used |21[]. In this algorithm, two jets are combined in to a single jet if 

2min(E t 2 ,Ef)(l-cos%) 
VDij = (59) 

is smaller than the jet resolution parameter yj. For pure QCD events, this algorithm 
tends to emphasize 2-jet events as compared to other algorithms and suited better for 



resummation purposes p3[ . When a photon in the final state is observed, we find higher 
1-jet rate and lower 2-jet rate and the QCD corrections are smaller as compared to the EO 
cluster algorithm. 



6 Conclusions 

Photon production in association with hadrons in e + e~ annihilation provides us interesting 
information on the non-perturbative component of the photon and new possibilities to test 
the underlying structure of perturbative QCD. 

In this paper we paid special attention to the importance of the correct treatment of 
the collinear photon-quark region. It was shown that next-to-leading and higher orders 
the perturbative part can only be defined using some non-physical parameter, no matter 
whether non-isolated or isolated photon production is considered. The physical cross 
section defined as the sum of the perturbative and non-perturbative part is, of course, 
independent of such a parameter. 

We briefly reviewed the theoretical description of the inclusive non-isolated photon 
production in e + e~ annihilation. It was pointed out that the LEP data can be used 
to constrain the parametrization of the fragmentation functions of the photon, D 7 / q (x, //), 
D 1)9 (x, fi). The measurement of these fragmentation functions would give important input 
information for the other inclusive photon production measurements at hadron colliders 
and at HERA. Furthermore, one could test the anomalous /^-dependence at asymptotically 
large /i values predicted by perturbative QCD. 

The case of isolated photon production was studied as well. Under well defined cir- 
cumstances, isolation can suppress the numerical contribution of the non-perturbative 
contributions. We pointed out that the non-perturbative ( "anomalous" ) contribution can 
be sizable only for E 7 > y / i/2/(l + e c ), where e c is the energy fraction in the isolation 
cone with respect to the photon energy. When a jet algorithm is used, then the non- 
perturbative contribution is expected to be further suppressed for isolated photon plus 
n-jet cross section for n > 1, but not for n — 1. 

We demonstrated the difficulty due to the quark photon collinear singularity with 
careful calculation of the next-to-leading order QCD corrections to isolated photon plus 
one or two jets. We argued that in the case of isolated photon plus 2-jet production indeed, 
as suggested by Kramer and Lampe, the perturbative contribution dominates the physical 
cross section. The next-to-leading order corrections are calculated by developing a Monte 
Carlo program which can be used to calculate the perturbative corrections to any physical 
quantity. 
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Appendix A 



(3) 

In this appendix we give the explicit expressions for the My and M„- expressions used 
in formulas fl34|) and ([!7]) respectively. We shall make the following renaming: 

q — > particle 1 

q — > particle 2 (60) 

g — > particle 3 

7 — > particle 4. 

Then the following relations are valid: 

M 23 = M 13 (l <-> 2), M 14 = M 13 (3^4), M 24 = M 13 (l ^ 2, 3 <-> 4), (61) 

therefore, it is sufficient to spell out Mi 3 . The corresponding expression can be obtained 
from the four-parton matrix element given in Appendix B of ref. Jl6| by setting Nq = 
and Tr = 0. After performing partial fractioning one obtains 



M 



13 



4:71 



2 



2yi22/i2s(l + Vu) 2yu(l - 2/24) 2(1 - 2/13)2/23 



,2/1342/234 (2/13 + 2/23) 2/234(2/13 + 2/24) 2/134(2/13 + 2/24) 
+ T2 - (2/242/34 + 2/122/34 + 2/132/24 ~ 2/142/23 + VliVls) + ' ' 1 



2/134 2/13 + 2/24 

2/122/242/34 + 2/122/142/34 " 2/132/24 + 2/132/142/24 + 22/122/142/24 / 1 1 



2/134 (2/13 + 2/14 + 2/23) V 2/13 + 2/23 2/13 + 2/14 , 

2/122/242/34 + 2/122/142/34 ~ 2/l 4 2/23 + 2/142/232/24 + 22/122/142/24 / 1 1 

2/234(2/13 + 2/23 + 2/24) V2/13 + 2/23 2/13 + 2/24, 

2/122/232/34 + 2/122/132/34 " 2/142/23 + 2/132/142/23 + 22/122/132/23 / 1 1 



+ 



2/134(2/13 + 2/14 + 2/24) V2/13 + 2/14 2/13 + 2/24 



22/122/1232/124 / 1 / 1 I 



2/13 + 2/23 + 2/14 + 2/24 V2/13 + 2/23 + 2/24 \ 2/l3 + 2/24 2/l3 + 2/23 

1/1 1 

+ : : : + 



2/13 + 2/14 + 2/24 V2/13 + 2/14 2/13 + 2/24, 
1/1 1 

+ ; ; ; + 



2/1342/234(2/13 + 2/24) 



2/13 + 2/14 + 2/23 \VJ3 + 2/14 2/13 + 2/23 / / 
(2/122/34 - 2/132/242/34 + 2/142/232/34 + 32/122/232/34 
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2/134(1/13 + Vu) 



+3J/122/141/34 + %i 2 ?/34 - 2/132/232/24 + 2j/i 2 2/232/24 
-2/132/142/24 - 22/122/132/24 + 22/i 2 2/24 + 2/142/23 
+ 22/122/23 + 2/142/23 + 42/122/142/23 + 42/i 2 2/23 

+22/122/14 + 22/122/132/14 + 4y 2 2 yi 4 + 2yj 2 y 13 + 2y? 2 ) 

(2/142/24 + 2J/142/23 + 22/122/14 + 2/ 2 3 
+2/132/23 + 22/132/24 + 2?/i22/l3 + Vu) 



Due to partial fractioning, this expression is finite if a single y^ —>■ (and for the same 
reason the expression is lengthy.) 

The virtual corrections can also be obtained easily from eq. (2.20) of ref. [TB[] by setting 
N c = and T R — 0. In our decomposition 



M (3) = 

99 91 



47T 5 



1 



1 

27 



(3 - 2 log 1/12) 



(62) 



and the finite part is 



M 



(3) 



1 



4vr 2 



2/12 



2/12 + 2/14 + 2/12 + 2/24 



2/12 2/12 + 2/24 2/12 + 2/14 



2/14 



2/24 



+ log 2/14 

+ log 2/24 

-2 



42/12 + 2 2/l22/l4 + 4J/122/24 + 2/142/24 



(2/12 + 2/24) 2 



42/?2 + 22/122/24 + 4i/i22/l4 + 2/142/24 



(2/12 + 2/14) 2 



2/l 2 2 + (2/12 + 2/14) 2 j-, / v . 2/l 2 2 + (2/12 + 2/24) 2 

— 2/24J H 



+ 



2/142/24 
2/l 2 4 + J/24 



2/142/24(2/14 + 2/24) 

2/24 2/14 2yi2 

2/14 2/24 2/142/24 7 V3 



2 log 2/12 



2/142/24 

2/i 2 2 



+ 



R(yi2,yu) 

22/12 



,(2/14 + 2/24) 2 2/14 + 2/24 



4 ! — + — — j ( ^7T 2 - log 2 J/12 - 



(63) 



where 



R(x,y) 



log x logy — log a; log (1 — y) — logy log (1 — x) H — 7r 2 

6 



Li 2 (x) - hi 2 (y) 



(64) 



and 



Li 2 (x) 



* log(l - g) 
dz 

z 



(65) 
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Figure captions 



Figure 1 Typical Feynman diagrams contributing to the calculation of the 0(a s ) correc- 
tions to the inclusive quark production in e + e~ annihilation. 

Figure 2 Leading order hard-scattering cross section for inclusive isolated photon produc- 
tion calculated in the MS subtraction scheme (eq. [27]) and in the cone subtraction 
scheme (eq. [H]). 5 C = 15° and e c = 0.1 isolation parameters were used. 

Figure 3 The dependence of the QCD corrections to the isolated photon plus 1-jet pro- 
duction on the unphysical parameter yo when physical cuts are applied (see formula 
(|44|)) — solid curves — and with unphysical cuts (only quarks are cut) — dashed 
curves. The slicing method in the pole approximation was used with y 7 = yj. 

Figure 4 The dependence of the QCD corrections to the isolated photon plus 2-jet pro- 
duction on the unphysical parameter y when physical cuts are applied (see formula 
(|44|)) — solid curves — and with unphysical cuts (only quarks are cut) — dashed 
curves. The slicing method in the pole approximation was used with y 1 = yj. 

Figure 5 Leading order photon energy spectrum of the partial widths T(Z — > 7 + njets), 
(n = 1, 2) normalized to lCT 3 r(Z — > hadrons) at yj = 0.1, 5 C = 15° and e c = 0.1. 

Figure 6 Next-to-leading order photon energy spectrum of the partial widths T(Z — > 
7 + n jets), (n = 1, 2) normalized to 10~ 3 r(Z — > hadrons) at yj = 0.1, 5 C = 15° and 
e c = 0.1. 

Figure 7 Partial widths T(Z — ► 7 + n jets), (n = 1,2) as a function of y,j normalized 
to 10~ 3 r(Z — > hadrons) when cone isolation of the photon is used (solid lines) with 
a s (Mz) = 0.12, a = 1/137. The dashed curves represent the scale dependence 
between scales \i = Mz/2 and fi = 2Mz- 

Figure 8 Partial width Y(Z — > 7 + 1 jet), as a function of yj normalized to 10~ 3 T(Z — > 
hadrons) (solid line) when invariant mass isolation is used with y 1 = yj, e c = 0.1, 
a s (M z ) = 0.12, a = 1/137. The dashed curves represent the scale dependence 
between scales /i = M^/2 and fi = 2Mz- The dashed dotted curve is the partial 
width calculated with y 7 = 0.005 and the long-dashed short-dashed curve is that 
with y y = 0.001. 

Figure 9 Partial width Y(Z — > 7 + 2 jets), as a function of yj normalized to 10~ 3 T(Z — > 
hadrons) (solid line) when invariant mass isolation is used with y 7 = yj, e c = 0.1 
,a s (Mz) = 0.12, a = 1/137. The dashed curves represent the scale dependence 
between scales \i = Mz/2 and \i = 2Mz- The long-dashed short-dashed curve is the 
partial width calculated with ?/ 7 = 0.001. 
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Figure 10 Distribution of the photon transverse momentum with respect to the thrust 
axis. The photon was isolated using invariant mass isolation with y 7 = 0.005 and 
0.06. The dotted histograms show the scale dependence between scales /x = M z /2 
and /x = 2M Z . The width of one bin is M 2 /100. 

Figure 11 Partial widths T(Z — > 7 + njets), (n = 1,2) as a function of y D normalized 
to 10~ 3 r(Z — > hadrons) when Durham clustering algorithm and cone isolation of 
the photon is used (solid lines) with a s (Mz) = 0.12, a = 1/137. The dashed curves 
represent the scale dependence between scales \i = Mz/2 and /x = 2Mz- 
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